Collective modes of monolayer, bilayer, and multilayer fermionic dipolar liquid 
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Motivated by recent experimental advances in creating polar molecular gases in the laboratory, 
we theoretically investigate the many body effects of two- dimensional dipolar systems with the 
anisotropic and 1/r 3 dipole- dipole interactions. We calculate collective modes of 2D dipolar systems, 
and also consider spatially separated bilayer and multilayer superlattice dipolar systems. We obtain 
the characteristic features of collective modes in quantum dipolar gases. We quantitatively compare 
the modes of these dipolar systems with the modes of the extensively studied usual two-dimensional 
electron systems, where the inter-particle interaction is Coulombic. 
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I. INTRODUCTION 

Recent experimental progress in producing and ma- 
nipulating ultracold polar molecules with a net electric 
dipole momentir— , provides new possibilities to explore 
novel quantum many-body physics in such systems^—. 
A series of theoretical work have been done within such 
dipolar system, such as the stable topological p x + ip y 
superfluid phase created with fermionic polar molecules 
with large dipole moment confined to a two dimensional 
geometry^, the existence of spontaneous interlayer su- 
perfluidity in bilayer systems of cold polar molecules^, 
the anisotropic Fermi liquid theory for the ultra cold 
fermionic polar molecules^, the zero sound mode in three 
dimensional (3D) dipolar Fermi gases^ 2 -, the superfluid 
properties of a dipolar Bosc-Einstcin condensate (BEC) 
in a fully three dimensional trap^ and the finite temper- 
ature compressibility of a fermionic dipolar gaa^. 

Most physical systems have long-lived excited states by 
conserving the total number of particles. These excita- 
tions have a bosonic character and are known as collective 
modes. A collection of charged particles is characterized 
by a collective mode associated with the self-sustaining 
in-phase density oscillations of all the particles due to the 
restoring force on the displaced particles, which arises 
from the self-consistent electric field generated by the lo- 
cal excess charges^. A two dimensional (2D) charged 
system can support density oscillations with the long- 
wavelength dispersion uj oc Jq. A similar collective 
mode occurs in a neutral Fermi system^. A repulsive 
short-range interparticle potential is sufficient to guaran- 
tee such a mode. This resulting density oscillation turns 
out to have a linear dispersion relation cj oc q in the long 
wavelength limit and is known as zero sound. The zero 
sound mode in connection with the RPA linear response 
theory was originally discovered by Landau as a collec- 
tive oscillation of the Fermi liquid with short-range inter- 
particle interaction as appropriate for a neutral system. 
Thus the zero sound and plasma oscillation are physically 
very similar, both being collective modes of an interact- 
ing Fermi liquid. However, the zero sound is physically 
very different from ordinary first sound, despite the sim- 



ilar dispersion relations u> oc q. 

In this paper we investigate the collective mode of the 
dipolar system, which has anisotropic and 1/r 3 dipolar 
interaction instead of the isotropic 1/r Coulomb inter- 
action. The collective mode of 3D dipolar systems is 
recently discussed as a solution of the linearized Boltz- 
mann equatio n 11 ' 12 . However, our theory is based on 
the leading order expansion of the dynamically screened 
dipolar interaction, the so-called infinite bubble diagram 
expansion, with each bubble being the noninteracting ir- 
reducible polarizability. In our approach, the detailed 
form of the long-wavelength dispersion is fixed by the 
behavior of q 2 V(q) as q — > 0, where V(q) is the dipolar in- 
teraction in momentum space. Since V(q) is anisotropic 
and behaves as a short-range potential we have the very 
interesting and zero sound like-collective modes in the 
long wavelength limit. We also consider a double layer 
dipolar system formed by two parallel single-layer dipo- 
lar systems separated by a distance a and a multilayer 
dipolar superlattice made of periodic arrays of 2D dipo- 
lar systems in the direction transverse to the 2D plane. 
Collective modes of 2D multi-layer structures have been 
extensively studied since the existence of an undamped 
acoustic plasmon mode was predicted in semiconductor 
double quantum well systems^. The collective modes 
can be detected with experimental probes that couple 
directly to the particle density operators. Typical exper- 
iments for solid state systems are inelastic light scatter- 
ing spectroscopy^— , frequency-domain far-infrared 2 ^ or 
microwave spectroscopy, or inelastic electron-scattering 
spectroscopy 2 ^—. 

The layout of the paper is as follows: In Sec. [Til we de- 
rive both analytical formula and numerical results for the 
plasmon dispersion relation in the 2D monolayer dipolar 
gas. In Sec. IIIH we study the plasmon modes in bilayer 
dipolar gas and their loss function (spectral strength). 
In Sec. HVi we present analytical results of the plasmon 
modes in the dipolar superlattice system within a simple 
model. In Appendix [X] and [B] we provide the detailed 
calculation for the interlayer dipolar interaction and two 
summations in Sec lIVi respectively. 
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II. COLLECTIVE MODE IN MONOLAYER 
DIPOLAR SYSTEM 

We start from the fundamental many-body formula 
defining the collective mode of a fermionic dipolar sys- 
tem. The collective mode of a fermionic dipolar system is 
given by the dynamical structure factor S(q,u)), which is 
proportional to Im[e(q, w) _1 ], where e(q, uj) is the dynam- 
ical dielectric function of the system. The longitudinal 
collective-mode dispersion can be calculated by looking 
for poles of the density correlation function, or equiva- 
lently, by looking for zeros of the dynamical dielectric 
function. 



e{q,u) = l-V(q)IL (q,oj), 



(1) 



where q and u are, respectively, the 2D wave vector par- 
allel to the plane and the frequency, V(q) is the inter- 
action between dipolar molecules in wave-vector space, 
and Ho(q, to) is the leading-order irreducible polarizabil- 
ity (i.e., the so-called bare bubble or the Lindhard func- 
tion in the relevant dimension). 

The interaction between the dipolar molecules is spa- 
tially anisotropic, which depends not only on the distance 
between two dipole molecules but also the angle between 
their relative vector and dipolar orientations. For the 
dipolar system in a 2D plane (xy-plane), the interaction 
between two dipoles located at r\ and r 2 , respectively, 
within the layer can be written as: 



V 2 D(ri - ri) = 



ri - r 2 \ 



r(l - 3 sin 2 6» B cos 2 0) (2) 



where d is the electric dipole moment. If the external 
electric field E is set in the xz-plane with an polar an- 
gle of 8 E , then <j> is the azimuthal angle relative to the 
x-axis. The configuration is depicted in Fig. [1] After 
Fourier transformation, we could get the dipolar interac- 
tion in the wave vector space. In order to handle with 
the short distance divergence of the 2D dipolar interac- 
tion and since we are more interested in the long wave 
length limit, we use the short distance cutoff c beyond 
which the dipolar interaction, given by Eq. [2] is valicUi. 
For the dilute Fermi gas, the short distance cut-off is 
set to satisfy the relation kpc <C 1. Then we have the 
interaction in wave vector space 

V 2D {q) = 2TTd 2 P 2 (cos9 E )(- - q) + nd 2 qsm 2 E cos2<f> q , 

c 

(3) 

where P 2 (cos8p) is the second Legendre polynomial. Al- 
ternatively, if we start from the 3D dipolar interaction 
and assume a fixed Gaussian density profile in the z direc- 
tion n(k z ) = e~ k * w / 4 , where w characterizes the typical 
confinement size of the two dimensional bipolar system 
in the z direction, then integrating over the z direction, 
as shown in the Eq. 3 of Ref. [Tj], yields the same re- 
sult with a numerical factor of order unity in front of 
l/w. From Eq. we can see that the first term of the 




FIG. 1: (Color online). Schematic sketch of the 2D monolayer 
dipolar system in momentum space. The polarization direc- 
tion of dipolar molecules is controlled by the external electric 
field E, which is in the xz plane. 



V 2 o(q) is isotropic, which can be either positive or neg- 
ative depending on the direction of the external electric 
field. While the second term of Eq. ^ is the anisotropic 
component, which can also be cither positive or negative 
but depending on the direction of q. 

The leading order irreducible polarizability is given by 
the bare bubble diagrams^ and for the spinless fermions 
we have 



n(q,w) = 



d 2 k n F (£ k ) - n F (£ k ^ 



(2tt) 2 



Ck - £k+q + If] 



(4) 



where £ k = h 2 k 2 /2m is the energy of single particle and 
rip is the Fermi distribution function. At zero temper- 
ature, the polarizability was calculated by Stern 2 ^ and 
the exact expressions, written in terms of dimensionless 
parameters x = q/2kp (hp is the Fermi wave vector) and 
8 = uj/AEp (Ep is the Fermi energy), is given by2£: 



RelRz, 8) = sgn(x - 



6 \2 



1 
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lmU(x,S) 



m 

47T.T 



V X 



5 \2 



i-(*-§) 2 



(5) 

where RcII(a;, 8) and ImII(a;, 8) are the real part and the 
imaginary part of the polarizability, respectively. 

We could further introduce some dimensionless param- 
eters which we will use in the following discussions: s = 
and r = 2d 2 kpm/H 2 which describes the strength of 
dipolar interaction. For example, experimental parame- 
ters of fermionic 40 K 87 Rb (to = 127 amu,d = 0.57 Dcbye, 
n = 10 8 cm~ 2 ), we have c = 10 nm (ck F ~ 0.03 <C 1), 
s ~ 31 (s>l) and r = 4.4. We use these parameters in 
the numerical calculation. 

For a single layer fermionic dipolar system, the collec- 
tive modes can be found by looking for the zeros of the 
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dynamical dielectric function, i.e., 

e(g,w) = l-V(g)II(g,w) = 0. 



(6) 



When the direction of electric field satisfies 3 cos 2 9 = 1 
(or 9 ~ 55°), the short distance cut-off disappears and 
the plasmon dispersion relation becomes 



First, we consider the leading order wave vector depen- 
dence of the collective mode. In the long wavelength limit 
(lu 3> qvF, where vf = kp/m is the Fermi velocity), the 
2D polarizability becomes 



n(g,w) 



q 



q 



a— + o{ij), 

ui z a; 4 



(7) 



where a = n/m, (n is the 2D density of the dipolar 
molecules and is related to the Fermi wave vector fcp 
as kp — %/ 47m). Then, we have the collective mode in 
the long wave length (q — > 0) 

uj(q) ~ q\l (27rd 2 P 2 (cos 9 E )(- - q) + nd 2 qsm 2 9 E cos 2<f> Q )- 

c \ 

(8) 

As 9e varies from to \ by changing the direction 
of external field, the dipolar interaction changes from an 
isotropic repulsive to an attractive one at large value of 
8 p. The plasmon mode given in Eq. depends on both 
the direction of the external field and the direction of 
the 2D wave vector. Since the 2D dipolar interaction 
in the wave vector space has both s-wave and d-wave 
symmetry as shown in Eq. <j3j> , we only need to consider 
the case with <j) q in the range [0, -|]. As special cases 
we investigate the plasmon modes along the x-axis and 
y-axis, which corresponds to <p q = and <j> q = j. For 
<fi q = 0, we have the plasmon mode at long wave length 
limit from Eq. (|8]l : 

q\l (2ird 2 P 2 {cos9p){- - q) + Trd 2 qsm 2 9 E )-. (9) 
V c m 

We find from the above formula that when the electric 
field is perpendicular to the xy plane (i.e. 9p = 0), the 
plasmon mode becomes 



u(q -4 0) ~ q 



2Trd 2 n 



Using kp = V 47T?i and s = 



md 2 
Wc 



we have 



){q -4 0) ~ VF\fsq. 



(10) 



(11) 



Since s>l, the plasmon mode shown in Eq. [10] satisfies 
the consistency criterion w 3> g«F, which has been used 
to get the 2D polarizability up to the leading order in 
wave vector. Note that for w > vpq the mode lies above 
the single particle excitation (SPE) regime (or particle- 
hole continuum) and prevents its direct decay through 
coupling to particle-hole continuum. The SPE region is 
defined by the nonzero value of the imaginary part of 
the total dielectric function, Im[e(g, ui)} ^ 0, which gives 
rise to the damping of a plasmon mode by emitting a 
particle-hole pair excitation 3 -!. 



<r 2 



3m 



rq 



-v F q. 



(12) 



We see that the undamped plasmon mode in Eq. (Tl2j) 
exists only for q > q c = Ylkpjr. For q < q c the mode 
enters into the single particle excitation region and it is 
damped by producing particlc-hole pair. As the external 
electric field is further tilted leading to 3 cos 2 9p < 1 the 
solution for u> satisfying Eq. ^ is purely imaginary, and 
there is no well-defined collective mode. Thus, along x 
axis, the direction of the external electric field must be 
smaller than the critical direction 9 C in order to exist an 
' undamped plasmon mode. 

Now we consider the collective mode along the y-axis 
(i.e. <f) q = 5). The plasmon mode for this case derived 
from Eq. ([8]) is given by 

U(q) ~ q\ (2nd 2 P 2 (cos 9 p)(- - q) - nd 2 q sin 2 9p) — 

V C 771 

(13) 

For an electric field being perpendicular to the single 
layer plane 9 E = 0, the plasmon mode of the system 
is isotropic and it is the same as given in Eq. (fTTj) . For 
3 cos 2 9 p = 1, there is an undamped mode at u> <C qvp 
and q > 2k f because Lindhard function becomes nega- 
tive and the dipolar interaction V2i)(g) is also negative 
along y-axis for 3 cos 2 9 E = 1. To get this unusual mode 
we expand the 2D polarizability for uj -C qvF as 



n(M) = f 



- 1 - 



(x 3 (-i + x 2 ) 3/2 ) 



(14) 

Using this large q behavior of the polarizability and 
Eq. we have the undamped low energy collective mode 
for 4> q = -| and 9 E = 6, which is consistent with the nu- 
merical result shown in the left panel of Fig. [5J 



v F q 
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The low energy mode at large wave vectors does not ex- 
ist in the usual two dimensional electron system (2DES) 
since the interaction of 2DES is isotropic and positive. 
This mode is unique for a fermionic dipolar system. Com- 
pared with the 2D fermionic dipolar system, the long 
wave-length plasma frequency for the extensively stud- 
ied 2DES is written as 3 ^ 



uj 2 (q -4 0) = 



— g 1 / 2 + 0( ? 3 ' 2 ) (16) 



where n and m are the charge carrier density and the 
effective mass of the charge carriers in the 2DES, respec- 
tively. The plasma frequency for the 2DES is isotropic 
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<f> q = Q 




the dielectric tensor is given by 
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FIG. 2: (Color online). Calculated plasmon mode dispersions 
of single layer dipolar system for different external electric 
field direction 6e- The shadowed region indicate the single- 
particle excitation (SPE) region. The interaction parameter 
r = 4.4 and s = 31 are used in this calculation. The mode 

TV 

dispersions (a) for <f> q = and (b) for cj> q — — are shown. 



and proportional to ^fq 1 while the plasma frequency of 
the 2D dipolar system is characterized by Eq. [8j which 
is anisotropic and has two different dispersions at long 
wave-length limit. 

In Fig. [21 we show numerically calculated collective 
mode dispersions of single layer dipolar system for three 
different external electric field directions Be- The cal- 
culated numerical results agree well with the analytical 
results discussed above. We choose Oe = 53° (it can be 
any angle close to the critical value 9) to investigate the 
behavior of the collective mode near to the critical angle 
Be = 9. 

When the 2D wave vector is along the x-axis (<j) q = 0), 
the slope of plasmon dispersion decreases as the electric 
field direction increases. In particular, for 9e > ~ 
55°, the plasmon mode enters into the SPE region and it 
is overdamped by producing particle-hole pair. On the 
other hand, when the 2D wave vector along the y-axis 
(<j> q = 77), the plasmon mode has similar feature as the 
case for <f> q = and for smaller angle of 9e- When Be 
approaches to the critical value 8 (i.e., 9 ~ 55°), the 
long- wave length plasmon mode approaches to the upper 
boundary of the SPE region (i.e., u> = q 2 /2m+VFq)- The 
main difference between (f> q = and <p q = tt/2 cases is 
that there is a short-wave length plasmon mode for the 
latter case even at 3 cos 2 9e < 1. This plasmon mode for 
4> q = n/2 appears below the lower boundary of the SPE 
region (i.e., u = q 2 /2m — VFq), where both the dipolar 
interaction and the Lindhard function are negative. 



III. PLASMON MODE IN BILAYER DIPOLAR 
SYSTEM 

For a bilayer system, which is parallel to the xy plane 
and is separated by a distance a, we need to consider 
the generalized dielectric tensor in order to find the 
collective modes. Within the RPA the Im component of 



ei m (q,u>) = Si m ~ Vi m {q)U m (q,ui), 



(17) 



where I, m = 1 or 2 with 1,2 denoting the index of two 
different layers. Then the plasmon modes are given by 
the zeros of the two-component determinantal equation, 



i.e. 



iq,u>) = [1 - Ku(?)IIi(ff,&;)][l - V 22 (q)U 2 (q,u)} 
-Vi 2 {q)V 21 {q)Ui{q,u)Tl 2 {q,u) = 0, 



(18) 



where 11; is the polarizability of layer I and is given by 
Eq. ([3]), and Vu and Vi m are, respectively, the intralayer 
and intcrlayer dipolar interaction. The intralayer dipolar 
interaction is the same as given in Eq. ((3)) 

Vu(q) = V 22 (q) - 2^ 2 P 2 (cos^)(i - q) 

9 2 C 

+7TGrsin 9ECOs2(j> q q , , 

1 ( i9 J 
= 7rd 2 (3cos 2 6 E -1)- 
c 

+2d 2 (77r(sin 9e cos 2 (j) q — cos 2 6e)- 

For the 2D bilayer dipolar system the intcrlayer inter- 
action V12 is given in real space 



V 12 (f) 



{r 2 +a 2 )i 



1 



3(r cos</>sin( 



h 2 



(20) 

where 9e is the polar angle of an external electric field 
which is applied in the xz plane and <p q is the azimuthal 
angle measuring from x-axis. Note that Vi 2 (—r) ^ Vi 2 (r) 
if the electric field is not perpendicular to the xy plane, 
which gives rise to the imaginary component for the in- 
tcrlayer interaction in the wave vector space. The de- 
tailed calculation for Vi 2 (q) is given in Appendix. [Aland 
we express the following explicit form of the interlayer 
interaction in momentum space as: 



Vu(q) = 2d 2 q-Ke- a i (sin 2 9 E cos 2 (j) q - cos 2 9 E 
— sm29E cos(j) q i) 



(21) 



and 



V 21 {q) = 2d 2 q-Ke- a i (sin 2 9 E cos 2 (j) q - cos 2 9 E 
+ sin29E cos <j) q i) 



(22) 



We note that both the interlayer and intralayer inter- 
action depends on both the direction of momentum q and 
the direction of an electric field E. In addition, there is 
a short distance cutoff in interlayer interaction, Eq. (|19[) , 
but not in the intralayer interaction. The interlayer dis- 
tance a should be larger than the short distance cut- 
off c beyond which the interaction between two dipole 
molecules can be described by the dipolar interaction. 
Typically in the cold atomic system, the interlayer dis- 
tance a — 500nm and we use c = lOnm as the short 
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distance cut-off for the intralayer interaction as done in 
the single layer dipolar system. Alternatively, we can as- 
sume that the two layers separated by a distance a are 
strongly confined in the z direction with Gaussian distri- 
bution in the z direction, n(z) oc e ^ (z±a ^ 2 ^ l 2w . In the 
wave vector space we have n(k z ) = e~ w k z±* d / 2k ~ . Then, 
we get the same interlayer interaction when we integrate 
out the k z part of the 3D dipolar interaction. The numer- 

a 2 

ical factor before l/w is suppressed by of Eq. (1) 

in Ref. [l4| , which means that the term with 1 /w can be 
neglected as long as a/w ^> 1. We numerically checked 
that the interlayer interaction with Gaussian distribution 
agree well with Eq. (f2"Tj) as long as a/c > 5, which is the 
case we consider in our work. 



A. Long wavelength plasmon mode in bilayer 
dipolar system 

In this subsection, we derive the analytical results of 
plasmon mode in the long wave length limit. We first 
consider the leading order wave vector dependence of the 
plasmon mode in the bilayer dipolar gas. At zero temper- 
ature, the two dimensional non-interacting polarizability 
has the following limiting forms in the high frequency 
regimes (i.e., w ^ qvp)r^'- 



plasma frequencies for the 2DES: 



2 4 

q q 
Ui(q,uj) = a.i—z + 0(—t) 

CO OJ^ 1 



(23) 



where i = 1 or 2 with 1,2 denoting the two different 
layers, and ai = nijm, ni is the dipolar molecule den- 
sity of the ith layer. Then, combining Eq. (|18I19[) and 
Eqs. (f2~T j) - (|23[) . we obtain the long- wavelength plasmon 
modes of the bilayer dipolar system: 



c*4 ~ y (aiVn + a 2 V 22 



(24) 



±y/ (aiVn — a 2 V 22 ) 2 + iaiaaV^^i; 



where cj+ (lu-) indicates the optical (acoustic) plasmon 
mode where the density fluctuations in each layer oscil- 
late in-phase (optical) and out-of-phase (acoustic) rela- 
tive to each other, respectively. Eq. [24] is also valid for 
the ordinary bilayer 2DES, while the interlayer and in- 
tralayer interaction is written as: 



2ne 2 {ni + n 2 ) 1/2 
V Km 



' 27re 2 n 1 n 2 
Km(ni + n 2 ) 



(26) 



While in the weak coupling (qa 3> 1) and long wave- 
length limit , the modes are simply the respective two- 
dimensional plasma frequencies of the two components^: 



27rn 1 e 



Km 



,1/2 



(27) 



27m 2 e 2 1 



/2 



Now we investigate the plasmon modes of this bilayer 
dipolar system in two regimes: the strong coupling limit 
(qa <C 1) and the weak coupling limit (qa S> 1). 



1. Strong coupling limit ga < 1 

When the external electric field is perpendicular to the 
bilayer plane, i.e., 9e — 0, the interlayer and intralayer 
interaction become for qa <C 1 



^ii( (? ) = ^ 22 ( (Z ) = 2^ 2 (i- (? ) 

c 



(28) 



V 12 (q) = V 2 i(q) 



-2TTd 2 qer aq 



Then the long wavelength plasmon modes become 



q 2 TTd 2 



Q l 



+ a 2 ± v (ai — a 2 ) 2 + 4q 2 c 2 aia 2 



(29) 

Because the short distance cut-off in the intralayer dipo- 
lar interaction is set to satisfy kp 1 2 <C - and n,2 = 
2md 2 kp 1 2 , it is clear to see that the above two plasmon 
modes also satisfy the self-consistent criterion u> qvp. 
These two modes lies outside of the SPE region, which 
are undamped and stable. Since the external electric field 
is perpendicular to the xy plane, the bilayer dipolar sys- 
tem is purely isotropic, which gives rise to the plasmon 
mode to be independent of the wave vector direction, <f> q . 

For 8p — w 55°, the system is purely anisotropic and 
the plasmon modes are given by 



V 11 (q) = V 22 (q) 



V 12 (q) = V 21 (q) = 



2ne 2 



27re 2 

q 



(25) 



which arc both isotropic. In the strong coupling (qa <C 1) 
and long wave-length limit, we have the following two 



9 d 2 nq 3 
u± - — 



(<X\ + a 2 ) cos 2(p q 



±Vtai - a 2 ) 2 cos 2 20 9 + 4aia 2 (2 + cos 2^) 5 



(30) 



When the 2D wave vector is along the 2- axis (i.e., 4> q = 
0), the lu- plasmon mode is overdamped since u;_ be- 
comes pure imaginary. But, the in-phase mode is well 



6 



3 2 



6 E = 0° 








a — 500nm 






-^SPE 



It} J 

■* 
3 2 



= 0° 




- n 2 = ni 




a — 50 nm 









i.O 0.2 0.4 0.6 0.8 1.0 

q/2k F 



i.O 0.2 0.4 0.6 0.8 1.0 

q/2k F 



FIG. 3: (Color online). Calculated plasmon mode dispersions 
of bilayer dipolar system for two different interlayer distances 
(a) a = 0.5/im and (b) a = 50 nm at 9e = 0°. We use n\ — 
7i2 = 10 8 cm~ 2 and r = 4.4 and s = 31 in this calculation. The 
shadowed region indicate the single-particle excitation (SPE) 
region. 



defined as long as r > 1 (in order to satisfy the self- 
consistent criterion oj 3> qvp) and the plasmon mode is 
given by 



(a>i + a 2 ) + \J (ai + a 2 ) 2 + 32ai« 2 



(31) 



However, along the y-axis (<p q = 7r/2), both w_ and w+ 
plasmon modes are pure imaginary, and there are no well- 
defined collective mode in the long wavelength limit. 



2. weak coupling limit qa 3> 1 

In weak coupling limit e~ aq ~ and the interlayer 
interaction between two layers is much weaker than the 
intralayer interaction, which give rise to the two inde- 
pendent plasmon modes of each layer. Thus we have 
uncoupled plasmon modes as 



uj 2 + = g 2 ai(2 7 rd 2 P 2 (cos6 l o)(i - q) + ird 2 sin 2 6 q) 
cj 2 = q 2 a 2 {2Trd 2 P 2 {cos6o)(- - q) + nd 2 sin 2 9 q). 



(32) 



The detailed discussion about the existence of plasmon 
modes in this weak coupling limit is the same as that for 
single layer dipolar gas. If the two layers have the same 
density and mass, = w 2 , i.e., two plasmon modes 
have the same dispersion relation. 



B. Numerical results of plasmon modes in bilayer 
dipolar system 

In this subsection, we show our numerical results of 
plasmon modes for the bilayer dipolar system with typ- 
ical parameters of 40 K 87 Rb. In Fig. [3] we show the cal- 
culated plasmon dispersions for an external electric field 
perpendicular to x-y plane {6e =0) with equal densities 
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FIG. 4: (Color online). Calculated plasmon mode dispersions 
of double layer dipolar system for two different interlayer dis- 
tance and (f> q at 6e = 53° with ni = n,2 = 10 8 cm -2 . The 
shadowed region indicate the single-particle excitation (SPE) 
region. The solid lines are solutions for e(q,u>) = for r = 4.4 
and s = 31. 



of ni = 112 = 10 8 cm~ 2 . When the separation of two lay- 
ers are much larger than the inverse of the Fermi wave 
vector (i.e. weak coupling limit), the interlayer interac- 
tion is negligible and the bilayer dipolar system behave 
like two independent single layers. Thus, two plasmon 
modes are almost degenerated and there is only one plas- 
mon mode showed up in the left panel of Fig|3J When 
the two layers are getting closer, the interlayer interaction 
become stronger. As a consequence, the density fluctua- 
tion in each layer is coupled and the degenerated plasmon 
mode is split into two plasmon modes, especially at large 
wave vectors, as shown in the right panel of FigJ3] 

In Fig. |?]we show the calculated plasmon dispersions 
for Be = 53°, (i.e., the direction of an external electric 
field is close to the critical angle of 6e = 55°). The 
plasmon dispersions are calculated with equal densities of 
m = ri2 = 10 8 cm -2 for two different interlayer distance 
(a = 500nm and 50nm) and two different directions of the 
wave vector (4> q = and tt/2) because of the anisotropic 
properties of the interlayer interaction. (Note that for 
8e = the plasmon modes are isotropic and independent 
of <f> q ). Since 9e is close to the critical angle the plasmon 
modes approach to the upper boundary of the electron- 
hole continuum. For <f) q = (i.e. the 2D wave vector 
is along x-axis) both in-phase and out of phase plasmon 
modes are located above the SPE region and they are 
undamped. As long as 9e > 0e we find two undamped 
modes even though the two modes are almost degenerate 
for weak coupling limit. When the layer separation is 
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FIG. 5: (Color online). Calculated plasmon mode dispersions 
of double layer dipolar system for two different Be at cf> q = 
with ni = 10 8 cm~ 2 (hp = V47™i), i2/^i = 0.5 and inter- 
layer distance a = 0.5/^m. The shadowed region indicate the 
single-particle excitation (SPE) region. The solid lines are so- 
lutions for e(q, to) = for r = 4.4 and s = 31. The left panel 
is for the electric field perpendicular to the xy plane Be — 0°, 
while the right panel is for the electric field titled at Be = 53°. 



small (i.e., strong coupling limits), the out-of-phase mode 
(lower energy mode) gets closer to the SPE region. For 
Be > Be the out-of-phase mode enters into the SPE and 
is overdamped. Thus, there is only one undamped mode 
for 9e > Be- For <f> q = ir/2 both modes have similar 
feathers as that for <f> q = but they are very close to the 
upper boundary of SPE region. When the 9e is greater 
than the critical angle we find no plasmon modes above 

q 2 

the upper boundary of SPE region (i.e., to = Vveq)- 

2m 

But there are plasmon modes appearing below the lower 
boundary of SPE region (i.e., ui = q 2 /2m — vpq) for <p q = 
tt/2 (the lower panel of Fig. HJ) . As the interlayer coupling 
gets stronger, the in-phase mode with higher frequency 
is pushed into the SPE region and only the out-of-phase 
mode with lower frequency survived. 

In Fig. we calculated plasmon mode dispersions of 
bilayer dipolar system with different densities (712/^1 = 
0.5) for two different 9e at cj) q = 0. Fig. [5ja) presents 
the plasmon modes for the isotropic interaction, i.e., the 
external electric field is perpendicular to the xy plane 
(6e = 0°). For a large interlayer distance (a = 0.5/im), 
the coupling between the two layers is very weak and they 
behave like two independent layers as shown in Eq. (f3"2|) . 
From Eq. (|3"I2]) . we know that, for Be = and weak 
coupling limit, both plasmon modes have linear disper- 
sion relation and the slope of the plasmon modes is pro- 
portional to their densities. Thus, two separated plas- 
mon modes are undamped. Fig. [SJb) shows the plasmon 
dispersion of the anisotropic bilayer dipolar system for 
Be = 53°, <j) q = and a = 0.5/im. There is only one un- 
damped mode w+ because w_ lies inside the SPE region. 

In Figs. |B]and [7] we show the calculated loss function 
(i.e., -Im[e(fcF, of the bilayer dipolar system for a 

fixed wave vector (q = Uf)- The loss function is related 
to the dynamical structure factor S(q,ui) by S(q,uj) oc- 
Im[e(fcj?,w) _1 ] and the dynamical structure factor gives 
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FIG. 6: (Color online). The density plots of the absolute value 
of the loss function -Imjeffcpjw) -1 ] of bilayer dipolar system 
with same densities in (cj,<fi q ) polar coordinate space for fixed 
density m = ri2 — 10 8 cm~ 2 , r = 4.4, s = 31, a — 500nm and 
9 E = 0°(left) and E = 53° (right). 



the direct measure of the spectral strength of the various 
elementary excitation a 31 ' 33 . We plot the loss functions in 
{io,<j) q ) space to describe the angular dependence of the 
plasmon modes. The loss function is an experimental ob- 
servable which can be measured with Raman-scattering 
spectroscopies. The plasmon modes exist when both 
real and imaginary part of the dielectric function equal 
zero (i.e., Rc[e(g,w)] = and lm[e(q, uj)} — 0), which 
leads the imaginary part of the inverse dielectric func- 
tion -Im[e(/sp, w) _1 ] to be a (5-function indicating an un- 
damped plasmon modes. On the other hand, the broad- 
ened peak in the loss function gives the damped plasmon 
modes. The damping in the plasmon mode is induced by 
particle-hole pair creation because we neglect any disor- 
der effects in this calculation. 

In Fig. HJa) we show the result when the applied elec- 
tric field is perpendicular to the plane. In this case both 
interlayer and intralayer interaction are isotropic. The 
solid circle gives the undamped plasmon modes, at which 
the loss function becomes singular. Since two plasmon 
modes are degenerate and they are independent of the 
wave vector direction <f> q , the plasmon modes appear as 
one solid circle in the (uj,<j> q ) space. For 6e = 53°, the 
two plasmon modes can also clearly be seen in Fig(6jb). 
For Be = 53°, both intralayer and interlayer interaction 
are anisotropic, so that the plasmon modes appear as el- 
lipses in the (oj,4> q ) space and they become much closer 
along the y-axis. 

In Fig. \7\ we present the density plots of the loss func- 
tion -Im[e(fcir, w) _1 ] of bilayer dipolar system for Be = B, 
where the short distance cut-off disappears from the in- 
tralayer interaction and the bilayer system becomes to- 
tally anisotropic. In Fig.[7]the plasmon modes are shown 
for interlayer distance (a) a = 0.5/^m (strong coupling 
limit) and (b) a = 50nm (weak coupling limit). There is 
no (5-function like peak in the loss function for q = kp in 
both cases, which correspond to damped plasmon modes. 
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where a is the supcrlatticc period (i.e., the separation be- 
tween adjacent layers in the z direction), sgn(Z — I') equals 
1(1 > I') or -1(1 < I'), and U v (q,w) = U(q,uj) is the ir- 
reducible polarizability of each 2D dipolar system given 
in Eq|4j The difference between intralayer and interlayer 
interaction in calculating the dipolar system is that the 
intralayer interaction in momentum space has the short 
distance cut-off c, while the interlayer interaction does 
not depending on this short distance cut-off. 



FIG. 7: (Color online). The density plots of the absolute 
value of the loss function -Im[e(fcF, i^) _1 ] of bilayer dipolar 
system with same densities in (u},(j) q ) polar coordinate space 
for fixed density m = na = 10 8 cm -2 , r = 4.4, s = 31, 8e = 
and h — 500nm(left) and a — 50nm(right). 



IV. PLASMON MODE IN DIPOLAR 
SUPERLATTICE 

In this section, we discuss the plasmon modes of a su- 
perlattice system which consist of the infinite number of 
parallel and equally separated two dimensional dipolar 
system. Due to the long range behavior of the dipolar 
interaction, it is necessary to couple all the layers, which 
changes the dielectric function in Eq. to an infinite 
matrix cquatioi 



,32. 



\5 ltV -Vi ll ,(q)lli,(q,u))\=0 



(33) 



where I (or l')= 0, ±1, ±2, are the indices of the layers 
and each layer is placed parallel to the xy plane with the 
position z = la in the z direction. Vi t i'(q) is the dipolar 
interaction between the Ith and the I'th layer, which can 
be written as 

1 = 1': 

Vu - V (q) = 7rd 2 (3cos 2 6 E -1)- 



Vw(q) 



-2d 2 g7r(sin 2 Q E cos 2 4> q — cos 2 9e)\ 



2d 2 qTre- a i\ l - l '\ (sin 2 6 E cos 2 <j> q 
—isgn(l — I') sin20 £ ; coscj) q ) 7 



cos 2 Be 



The plasmon modes of the infinite periodic system can 
be calculated from the self-consistent field method de- 
scribed in Rcf. 34 1 . Since the supcrlattice is perfectly 



periodic in the z direction with periodicity a, we assume 
the following ansatz m(q,uj) = no(q,ui)e zkzla , where a is 
the interlayer distance and I is the layer index. The k z 
introduced in the above ansatz labels the induced den- 
sity fluctuation in the infinite periodic system, which is 
restricted in the first Brillouin zone of the superlattice, 
i.e., < k z < — . Then the plasmon modes for the dipo- 
lar superlattice system are given by the roots of following 
equation 



(34) 



l-U(q,u>)V (q)-U(q,cj) ]T V w ( q ) e - ik ^- l > = 

(35) 

where Vo(q) is the intralayer interaction in momentum 
space while Vi-i'(q) is the interlayer interaction. With 
the help of formula given in Appendix. [B] Eq. ([33]) can 
be explicitly written as: 



1 - H(q,Lu)V (q) - U(q, uj)2d 2 qn [(sin 2 E cos 2 (j> q - cos 2 9 E )Si(q, k z ) - sin 29 E cos 4> q S 2 (q, k z )] = 0, 



(36) 



r 



where 



and 



Si(q, k z 



cos k z a — exp (—aq) 
coshaq — cos k z a 



(37) 



We investigate the long wavelength plasmon modes for 
the strong coupling qa ^> 1 and the weak coupling case 
qa <C 1 analytically in the following two subsections. 



(38) 
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A. Strong coupling case qa <C 1 



rewrite Eq. (|36|) for qa <C 1 as 



First, we consider the k z ^ case. With the asymp- 

2 

totic form of the polarizability Tl(q, uj) ~ we can 



1 - 27rd 



2 n 9 



2 r 



P 2 (cos^) 



+ (sin 2 9e cos 2 



cos Ve, 



aq 



1 — cos k z 



qsin29E cos( 



sin fc,a 



1 — cos k z 



= 0. 



Then we have the plasmon modes for the supcrlattice dipolar system in the long length limit 



2 2 ' 

Ul ~ q 



2mrd 2 \P 2 (cosB E ) 



aq 



+ (sin Be cos &„ — cos Be) — 

1 — cos fc 2 a 

I 



(39) 



(40) 



The plasmon mode in superlattice system is strongly de- 
pendent on the direction of the external electric field Be 
and the wave vector direction <p q as well as the wave 
vector k z that labels the density fluctuation. When the 
electric field is perpendicular to the plane, we find that 
the plasmon mode has linear dispersion relation to oc q. 
For Be = B, in which the interactions are anisotropic, the 
plasmon mode dispersion becomes 



3 / 2nird 2 , aq cos 2<p q 
^ V 3m 1 — cos k z t 



2V2 



cos 4> q sin k z a N 
f — cos k z a ' 



(41) 

In order for this mode to lie above the SPE region, i.e. 
uj > qvF, it is required to be r > 1. 

For k z = 0, the long wavelength plasmon mode can be 
calculated from Eq. to be 

q 2 2nixd 2 [~P 2 (cos Be) ^_ 2 (sin 2 Be cos 2 4> q — cos 2 Be) 



(42) 

In particular, the plasmon dispersion for Be = becomes 
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2nird 2 {a- 2c) 



(43) 



From Eq. (|43[) it is clear that only for a > 2c there is 
well defined plasmon mode with linear dispersion rela- 
tion. For the total anisotropic case (Be = Be), the plas- 
mon mode also has linear dispersion relation which can 
written as: 



Anird 2 cos 2cj) q 
3ma 



(44) 



Thus the plasmon mode only exists for certain angles 
such that cos 2cj> q > 0. 

Compared to the dipolar superlattice system, the 
plasma frequency for the superlattice system with 
Coulomb interaction can be written as^ 2 -: 



w(q -)• 0; k z = 0) 



47re 2 n\i/2 



(45) 



which has the precise character of the corresponding 
three dimensional plasmon (with a finite gap) in the long 
wavelength limit. 



B. Weak Coupling case qa ^> 1 
For the weak coupling (qa ^ 1) Eq. (f3l)|) becomes 

i-n(?,w)v (g) = (46) 

for all values of k z . Thus the long wavelength plasmon 
mode of superlattice system becomes the mode of the 
single layer dipolar system as discussed in Sec[II] Each 
layer has its own plasmon mode in the weak coupling 
limits as expected. 

V. SUMMARY AND CONCLUSIONS 

In summary, we have derived the collective plasmon 
modes in 2D dipolar Fermi liquids, and also considered 
spatially separated bilayer and superlattice dipolar sys- 
tems within the random phase approximation, which is 
valid for weakly interacting system. We have also calcu- 
lated the loss function for bilayer fermionic dipolar gas, 
which could be measured in experiments. The 2D dipo- 
lar interaction with both s and d-wave components in 
the wave vector space leads to several unexpected fea- 
tures in the plasmon modes such as undamped mode if 
dipole along the z-axis and the anisotropic plasmon dis- 
persion relation if the dipole along other direction than 
the z-axis. Our predicted plasmon modes is clearly dis- 
tinguished from the extensively studied two dimensional 
electron system. 

Future work ought to include higher order corrections 
to the polarizability and the finite temperature correc- 
tions, yet this must await the successful fabrication of 
quasi two dimensional dilute fermionic dipolar gases. We 
note that our theory for the collective mode spectra of 
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2D dipolar Fermi liquids should have considerable rele- 
vance for the recently made polar molecular gases Ref. 0- 
[§]. At low enough temperatures, T <C Tp where Tp is 
the Fermi temperature of the dipolar system, our theory 
should in principle describe the laboratory polar molec- 
ular fermionic systems, and the excitation spectra of the 
dipolar system should have clear signatures of the collec- 
tive mode spectra presented in our work. 
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= J ~ It " r2a : 0S ^ e ° exp W . r)drc 



(a 2 + r 2 )^ 



(A7) 



2d 2 que aq sin 29q cos <j> q i 



which is pure imaginary number. Finally we have the 
intcrlaycr interaction in momentum space 
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V 12 {q) = 2d 2 qiTe a9 (sin 2 E cos 2 (j) q -co& 2 E sm29 E cos(/> q 

(A8) 

We also find V2i(q) from V 2 i(r) 



Appendix A 

Below we provide the detailed calculation of the inter- 
layer dipolar interaction in the wave vector space. In the 
real space the interlayer dipolar interaction becomes 



V 12 (r) = 



(r 2 + a 2 )"- 



^ 3(r cos</>sin0E + hcos9 E )' 
r 2 + h 2 



(Al) 

where <f> q is the angle between momentum q and the x- 
axis. We can devide above equation into three different 
partss depending on the symmetry 



V 21 {r) 



(r 2 + a 2 y. 



3(r cos^sin^s — h cos Be) 



P 2 + yp, 



(A9) 

the interlayer dipolar interaction V 2 i{q) in the wave vec- 
tor space is given the complex conjugate of V\ 2 (q), which 
can written as: 



V 21 {q) = 2d 2 q-ne ° 9 (sin 2 E cos 2 9 -cos 2 8 E +sm29 E cos(/> q 

(A10) 



Vl 2 a) (r) 



[r* + a "2 



"(I 



3a 2 cos 2 8 
r 2 + h 2 



), (A2) 



Appendix B 



and 



V} 2 \r) 



(r 2 + a 2 )' 



-(3r 2 cos 2 </>sin 2 9 



(A3) 



-(6r cos^sin^oacos^o), (A4) 



(r 2 + a 2 )* 

where V 12 (q) = V$ (q) - V$ (q) - V$ (q) . Then we have 



^ g a cos t 

Vi 2 ] (v) = d 2 /o°° XT , n a2 tsl ex P(*9 • r)rdrd4> 



{a 2 + r 2 Y- 

„ , sin 2 On „ 

2qnd 2 e- ac i{ - cos 2 9 , 

aq 



(A5) 



where Jo(x) is the Bessel function of the first kind and 
q ■ f = qr cos(0 — <p q ). 



r(b),^ ja r oc r 27r 3r 2 COS^ </> sin' a 



(a 2 + r 2 ) 2 



exp(i<f ■ r)rdrd(j) 

(A6) 



1 



2d 2 qire~ aq sin 2 O ( cos 2 0„), 

aq 



Below we provide detailed calculation of two infinite 
sums in Sec. IIVI 



exp[— |Z|ag — ik z la] 
e-xjp(—laq — ik z la) 



Si(q,k z ) 



and 



S 2 {q,k z ) 



2Re[- 

1 — cxp(— Laq — ik z la) 
cos k z a — exp {—aq) 
cosh aq — cos k z a ' 



i sgn(7)exp[— \l\aq + ik z l 
exp(—laq + ik z la) 



(Bl) 



a 



2Im[ 



1 — cxp(— laq + ik z la) 
sin k z a 



cosh aq — cos k z a 



(B2) 
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